# Required functions: `srs.1prop`

weather=read.table('http://www.webpages.uidaho.edu/%7Erenaes/Data/temps2.txt',header=T,sep='')
set.seed(44676)
# install.packages("sampling") if you have not yet done it
library(sampling)
N=88; n=25
s=srswor(n,N)
sobs=(1:N)[s==1]
weather.srs=weather[sobs,]
View(weather.srs)
marchwarm=subset(weather.srs,weather.srs$MT>32)
marchcount=length(marchwarm)
   
# CI with x=15, n=182
srs.1prop('weather.srs',x=marchcount,n=n,N=N)